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Abstract 

We propose a box and ball system with a periodic boundary condition 
(pBBS). The time evolution rule of the pBBS is represented as a Boolean 
recurrence formula, an inverse ultradiscretization of which is shown to 
be equivalent with the algorithm of the calculus for the 2JVth root. The 
relations to the pBBS of the combinatorial R matrix of U' q (A$) are also 
discussed. 

1 Preface 

In many physical phenomena, physical quantities and time-space variables are 
continuous. Discretized models, however, are used sometimes for simplifica- 
tion and/or speedup of computation. In discretizing process, it is preferable 
to keep the mathematical structure (symmetry, conserved quantities, etc.) of 
the original continuous systems. For integrable systems, i.e. the systems which 
exhibit solitonic natures and are described by integrable nonlinear partial dif- 
ferential equations, several effective methods of discretization have been estab- 
lished and many discrete integrable systems have been proposed [[[]. An effec- 
tive method to create a cellular automaton (CA) with solitonic natures is the 
ultradiscretization Q which is a limiting procedure to discretize the dependent 
variables of a partial difference equation. 

In this paper, we consider an integrable cellular automaton called box and 
ball system (BBS)|| [|. The BBS has been denned as the time evolution of 
a finite number of balls moving in one dimensional array of infinite number of 
boxes. We extend the time evolution rule of the BBS so that the BBS can be 
defined when we impose a periodic boundary condition. We show that the time 
evolution rule of the periodic BBS (pBBS) is given by a Boolean formula, which 
turns out to be an algorithm to calculate the 2Ath root of a given number 
through inverse ultradiscretization. It is well known that many of the discrete 
integrable systems are equivalent to good algorithms such as the diagonalization 
of matrix (QR algorithm, etc.), convergence acceleration algorithm (e algorithm, 
etc.), Karmarkar algorithm, BCH-Goppa decoding. Our Boolean formula is 
an example of such correspondence between discrete integrable systems and 
algorithms. 

In section 2, we introduce pBBS and explain its time evolution rule. In sec- 
tion 3, we obtain recurrence formulae on Boolean algebra for the time evolution 
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rule of the simplest BBS with box capacity one. The correspondence of the 
Boolean formula with the algorithm for the computation of the 2iVth root is 
shown. In section 4, we establish a relation of the pBBS to the combinatorial 
R matrix model and extend the time evolution rule to the general pBBSs with 
various box capacities. Section 5 is devoted to the concluding remarks. 



2 Box and ball system 
2.1 Infinite BBS 

The BBS is a reinterpretation and extension of the filter-type CA proposed by 
Takahashi and Satsuma]^, 

Let us consider a one dimensional array of infinite number of boxes. At the 
initial time t = 0, all but finite number of the boxes are empty, and each of the 
rest boxes contains one ball such as: 



The time evolution rule of this system from time t to t+ 1 is given as follows. 

1. Move every ball only once. 

2. Move the leftmost ball to its nearest right empty box. 

3. Move the leftmost ball of the rest ball to its nearest right empty box. 

4. Repeat the above procedure until all the balls are moved. 

An example of the time evolution is shown in Fig. |l|. In the example, we see 
the solitonic behavior of balls. This is a general behavior of the BBS and, 
starting from an arbitrary initial state, we can prove that the state asymptot- 
ically evolves into the state that consists of only freely moving solitons. It is 
also known that the BBS has infinite number of conserved quantities^. By 
replacing an empty box by " 0" and a filled box by " 1" , we regard the BBS as a 
dynamical system of "01" sequence, namely, a CA. The relation between BBS 




Figure 1: An example of the time evolution of BBS. 
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and (classical) integrable systems, i.e. integrable partial differential and/or dif- 
ference equations, was established by the notion of ultradiscretization (UD)^j. 
In the approach of UD, the BBS is constructed through a limit of the partial 
difference equation, which is obtained by reduction of the discrete KP equation 
(Hirota-Miwa equation). Since the discrete KP equation is equivalent to the 
generating formula of the KP hierarchy, the BBS naturally inherits its integra- 
bility and shows solitonic natures. In this sense, the BBS is called an integrable 
CA. 



2.2 Periodic BBS (pBBS) 

In this subsection, we extend the original BBS to that with a periodic boundary 
condition. Let the BBS consist of N boxes. To impose a periodic boundary 
condition, we assume that the Nth box is the adjacent box to the first box. 
(We may imagine that the boxes are arranged in a circle.) We also assume 
that the box capacity is one for all the boxes. Since the evolution rule of the 
BBS requires the definition of the leftmost ball in the BBS, wc cannot apply the 
original evolution rule directly to a periodic system. Instead, we consider the 
following time evolution rule: 

1. Move all the balls to their conterminous right boxes if the boxes are empty. 

2. Forget the boxes to which and from which the balls were moved in the 
first step, and move all of the rest balls to their conterminous right boxes 
if they are empty. 

3. Repeat the above procedure until all the balls are moved. 

An example of the movement of the balls by this rule is illustrated in Fig. ^. 

We see that this rule is equivalent to the original time evolution rule when 
we apply it to the BBS with an infinite array of the boxes. In this rule, the 
number of balls which move at each stage does not change in time evolution. 
We will denote by a periodic BBS (pBBS) the BBS with a periodic boundary 
condition which evolves in this rule. 
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Figure 2: An example of the modified rule for pBBS. 



3 



Example 1 (pBBS (N = 7)) 

Time evolution patterns of pBBS are shown in Fig. ^. (a) The initial state is 
" 1110000" . The time evolution pattern has the fundamental cycle 7. (b) The 
initial state "1101000" has the fundamental cycle 21. 

The time evolution of pBBS is represented by a diagram with N nodes and 
M lines each of which connects two of the distinct nodes as shown in Fig. ||. 

Each node corresponds to a box of the pBBS and a line is drawn from the 
box which contains a ball to the box to which the ball is moved. From the time 
evolution rule of the pBBS, we find the following properties of the diagram: 

• There is no intersection of the lines. 

• There is no line over the nodes which are not connected by lines. 

• If we remove the unconnected nodes, the diagram is divided into a set of 
blocks with even number of nodes which are connected by lines pairwisely. 

For convenience of explanation in the subsequent sections, we define indices 
{7t( n )}^=i °f t ne nodes in the diagram at time step t: 

• If a ball in the nth box is moved to the n + kth (modulo N) box at t + 1, 

jt(n) = k(e Z> ). 

• If a ball is moved to the nth box from the n — fcth (modulo N) box at 
t + 1, 7t(n) = -k. 

• Otherwise, 7*(n) = -co 

In the example shown in Fig. ^, the indices {7*(n)} = { — 1, —3, — oo, — oo, 9, 5, 
1, —1, 1, —1, —5, 1, —1, —9, — oo, 3, 1 }. The following Proposition is a direct 
consequence of the properties of the diagram. 
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(a) time (b) 

Figure 3: Examples of pBBS (N = 7). 
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Figure 4: A diagram of time evolution and corresponding time evolution of the 
pBBS. The integers below the nodes are the indices 7t(n). 

Proposition 1 

Let k and I ( I < k ) be positive integers. Then 



(a) If 7t(n) = k, then 7 t (n + k) 
1 < j < k - 1. 



-k, and -(k - 2) < 7t(n + j) < k - 2 for 



(b) Furthermore if there exists j such that 1 < j < k — 1 and 7t(n + j) = —I, 
then I < j and 7t(n + j — I) = I. 

In Proposition |], the indices are understood with the convention of modulo 
N, i.e. 7t(n) = 7t(n + iV). Hereafter, we often use this convention. 

There are several equivalent evolution rules for the pBBS. For example, 

1. At each filled box, create a copy of the ball. 

2. Move all the copies once according to the following rule. 

3. Choose one of the copies and move it to its nearest right empty box. 

4. Choose one of the rest copies and move it to its nearest right empty box. 

5. Repeat the above procedure until all the copies are moved. 

6. Delete all the original balls. 

It is not difficult to prove that the result does not depend on the choice of the 
copies at each stage, and that this rule gives the same time evolution pattern 
of the previous rule. An advantage of this rule is that it is straightforward to 
extend the rule to the BBS with many kinds of balls and various box capacities. 
We can also use the combinatorial R matrix of U' q (A^'_ 1 ) to give an equivalent 
evolution rule, the detail of which is presented in section 4. 



3 Recurrence equations and corresponding al- 
gorithm 

3.1 Boolean formulae of pBBS 

We show that pBBS introduced above can be formulated by Boolean algebra. 
Let N be the number of the boxes. The space of the states of pBBS is naturally 



5 



regarded as ¥ 2 N . We denote a state X(t) of the pBBS at time t by X(t) = 
(xi(t),x 2 {t), . . . ,XN{t)) € ¥ 2 N , where Xi(t) = if the ith. box is empty and 
Xi{t) = 1 if it is filled. Let A, V, ©, be AND, OR and XOR respectively. These 
Boolean operators are realized in ¥ 2 N as the map: ¥ 2 N x ¥ 2 N — > ¥ 2 N . For 
X = (xi,X2, ■ ■ ■ ,xn), Y = (3/1,2/2, ■ • ■ ,3/jv), they are defined as 



(X A Y)i 
(X V Y)t 

(x © Y)i 



XiVyi = XiVi + x l + y l 
Xi ffi 2/i = X{ + 3/j . 



We also define rotate shift to the right S: 

SX = (xn,xi,x 2 , ■ ■ ■ ,xn-i)- 
Next theorem gives an expression of T 

T : X(t)^ X(t + 1) 
in terms of these Boolean operators. 
Theorem 1 

Suppose that X(t) e ¥ 2 N is the state ofpBBS at time step t. We consider the 
following recurrence equations: 

A™ =X{t), =SX{t), (1) 



B (n+1) ._ g(A(n) A _g(n)) 



(n = 0,1,2,...). (2) 



Then, 



l(Hl) = A (A,, ®I(t), andB^^O, (3) 
where := (0,0,... ,0). 



Proof 1 

We define T>« = (df , df , ... ,d$) £ ¥ 2 N (i = 1, 2, . . . , N) as follows. If the 
Ith box contains a ball at time step t and it moves to the l + ith box att+l, then 
d[ l) = 1, otherwise df ] = 0. If we use the indices j t (l) defined in the previous 
section, the dehnition of £)W j s rewritten as 

d^f 1 — 1 if and only if 7t(Z) = i. 

Clearly, these give the unique decomposition of the initial state X(t) 
and the final state X(t + 1) as 

N 

X(t) = 0I>W, =0 (i^j), 

»=i 

AT 

TX(t) = S* l T> w , S i D^ A S^'dW = (i 7^ j), 

i=l 
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where 0^ D« := L> (1) ® ... D^. 

To prove the theorem, it suffices to prove the following formulae 



A< n > = !Wffi05 l D (,) , (4) 
i=l 

73 

£-„-!£(„) = x(i)g0flW, ( 5 ) 

i=l 

Indeed, if (Q) and (Q) noid, we Lave 

AW©x(t) = x(t)e0s*D (i) ex(t) 

»=i 

AT 
8=1 

= *(t + i), 

= (X(t) © X(t)) 

= 0. 

We prove and by induction. 
For n = 1 , 

= X(t)VSX(t) 



N 



= X(f)VS5 (1) V0SD (i) . 

Suppose that 0^ 2 SD®VX(t) ^ X(t), there exists £>(«) (q > 2) such that its 
n' component d^) = 1 and 7t(n' + 1) < — 1. However, hy the definition of D^ q \ 

a$ = 1 implies 7t(n') = g(> 2). From Proposition [1| (a), — g + 2 < 7t(n' + 1) < 
q — 2 and we find — q + 2 < 7t(n' + 1) < — 1, which contradicts Proposition |I] 
(b). Hence 0^ 2 SD® V = X(t). Since SD^ A X(t) = foiiows from 
TX(t) A X(t) = 0, wc Gnd 

A {1) =X{t)®SD {1) . 
Similarly, from the relations X{t) A 0^ 2 SD® = 0^ 2 SD^ and SD^ A 
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X(t) = 0, we have 

5-2^(1) = g-i A SX(t)) 



S- 1 ^X(t) A05I)wj 
S^ 1 (x(t) A 5L> (1) ) 8 S -1 (x(t) A SD ( 
(S^O)©^ 1 (&SD<A 



N 





N 

= 0^ } 

i=2 

Hence (Qj and (g) ioJci for n = 1 . 

Assume that (||) arid |5|) is true for n — k, then 

A {k+1) = ^X(t) © ^D (i) ^J V ^ fc+1 X(<) © S k+1 D^>j 
= ( X (t) © SWiA V ( S fc+1 £« j . 

\ i=l / \j=fc+l / 

By the definition of D^ k+1 \ we have 

(fc \ fc+1 

X(t) © S l D^ V 5' fe + 1 Z3(' £ + 1 ) = X{t) © S*I>W. 

Suppose that 

(k \ N k 

X(i)ffi0S* l £>« V 5 fc+1 £> w ^X(t)©0S ,l L>W. 
i=l / i=fc+2 i=l 

Noticing the fact that, for X(t) © 0* =1 = (af°, o%° \ ... , 0$), af ) = 1 

if and oniy if — fc < 7t(j), we find that there is at ieast one (g > fc + 2) one 
of whose components satisfies d^) = 1 and 7t(n' + fc + 1) < — fc— 1. However, by 
the definition of D^ q \ dffi = 1 implies 7t(n') = q(> fc + 2). From Proposition [I] 
(a), -q + 2 < 7 f (n' + fc + l) < q-2 and we find ~q + 2 < -f t (n' + k + l) < -k-1, 
which contradicts Proposition^ (h). Thus, we have 

fc+i 

A (fc+i) =x(t)©0S ,i £)( i ). 

i=l 

In a similar manner, we also obtain 

fc+i 

s -k-2 B (k+l) = X ^ e Q 
i=l 



Hence, from the assumption of induction, the formulae (Q) and (||) are proved 
to hold for < n < N — 1, which completes the proof of the thorem. 

This recurrence equation (j^) is expressed with only three operations, AND, 
OR and SHIFT, and has a simple form. The SHIFT operator introduces the 
right-and-left symmetry breaking that comes from the definition of the direction 
of the movement of balls. 

Since we immediately obtain the following Corollary: 

Corollary 1 

Suppose that X{t) £ ¥2 N is given as the state at time t. Then the state at next 
time X(t + 1) = TX(t) is calculated by recurrence equation as follows. 



A® := X(t), B (0) := SX(t) 
B (n+1) ._ S 2( A (_n) A _B(n)j 



(6) 
(7) 



and 



x(t + 1) = a^ n ^ e x(t), b(W2J) = o. 



(8) 



3.2 pBBS and numerical algorithm 

The formulae for time evolution of the pBBS (||) have simple and symmetric 
form, and we expect that they have some relations to a good algorithm. In this 
subsection, we show that they have indeed the same structure as that of the 
algorithm to compute Nth root of a given number. Henceforth, let the truth 
value "O(false)" and "l(true)" be equivalent to the integer 6 Z and 1 e Z. 
Then we can replace A and V with min and max as 

x Ay <==> min [x, y] 
x V y max [x, y] 

Following the notation in the previous section, we define that max and min 
act on Z N bitwisely. Then, Eq. (|J) can be rewritten by the equation of integers 
as 



^(n+i) _ max [A( n \B^] 
B ( n +X) = Smin [aW,BW] 



(9) 



We construct the difference equations corresponding to (|^) by means of inverse 
ultradiscretizationQ. Noticing the identity: 

max[x,y}= lim elog ( e x/t + e y/€ ) (x, y 

e— >+0 V / 

and min [x, y] = — max [—x, —y], we think of the difference equations: 
4" +,, = {a<"> + *,<'''} /2 

^-timr+wr 1 '- 1 {1£,£N> - (10) 
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The relation between @ and (10) is obvious. When we replace and b[ n ^ 
with e^ A( )S)i t e and e*- 5 ' respectively, and take a limit e — > +0, we obtain 



(|9|) from fllpp . The factor 2 in (|10|) is so chosen that the recurrence formulae do 
not diverge at n — > oo. 

When we disregard the space coordinates i in ([l0|), or consider the case 
N = 1, we have recurrence formulae 



a 



(n+1) _ fl(") + 6 (?t) 



a {n) + b {n) 

which is the well-known arithmetic-harmonic mean algorithm and we have 
lim a (n) = lim = Va^M ). 



The recurrence formulae ( |10[) for general AT is also considered as a numerical 
algorithm to calculate the 2iVth root of a given number. To see this, first we 
note that (^) has a conserved quantity C with respect to the step n, 



N N 



CW := [] <4 n) &f 5 - C {n ^ = • ■ • = = [] <4° ] bf ] ee C, (12) 
i=i i=i 

where |a| ', are the initial values. Then we can show the following Propo- 
sition: 

Proposition 2 

If all the initial values {a[°\ fy° } are positive, then they converge to the same 
value 



lim ai^ = lim fo£ = 2 x 

n — >oo n — >oo 



A' 



v U a ^ b{ " )= 2 ^ (for all k). 
\ i=i 



Hence, the recurrence formula of pBBS is regarded as a numerical algorithm of 
2Nth root. 

To prove the Proposition, we need a Lemma: 
Lemma 1 

For m > 0, a > 0, and e > which satisfy 

(2a 2 + 5a + 4)e < m, (13) 

if it holds that m — e < a[ n+1 ^ < m + ae, m — e < a\ n ^ and m — e < b[ n \ then 
a\ n \ Q < m + (2a + 2)e. Similarly, if it holds that m — e < frj+i < m + cue, 
m — e < a\ n ^ and m — e < b\ n \ then a\ n \ b[ n ' < m + (2a + 2)e. 



This Lemma is proved from ( 
give the proof of the Proposition 2| 



lOj) by straightforward calculations. Now we 
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Proof 2 

Let mW 



•4=1,. 



,N 



we have from ( |Ig| ) 



a„ (n) A (n) 



an 



d AfW := 



maxi=j. 



,(») l(«) 



Since 



(n) , (n) 



(n) ,(n) 



< a. 



(n+l) 



< 6 



(n+l) 
i+1 



< max 

< max 



,(«) /.(") 

4 i ' u i 



(n) , (n) 



we obtain 



< TO («+ 1 ) < M < M W for n = 0, 1, 2, . 



(14) 



From (R), we Bnd^m, m = lim„^oo m/™) and 3 M, M = limn—^ M^ n \ Clearly 
m > mW, M < M(") and to < M. 

We win" prove m = M by leading contradiction to the assumption to < M. 
Since and converge to to and M respectively, for all e > 0, there 

exists no such that < m — m/™' < e and < — M < e for v n > n . 

We take £ = min [(M - m)/2 N+1 , m/2 2N+i ] . Note that, with this choice, the 



(n +N) 



inequality ( |Ti[ ) holds for < a < 2 N+1 . For n = no + N , there exists a-j ' or 

rnS na+N > , noticing vn — e < 



bf 0+N) which is equal to m n °+ N . When a ( f 0+N) 

m (n +N) < m, m — e < a\ n ° +N 1 ' and to — e < b\ n ° +N l \ we obtain from 

(no+W-l) ,(n +JV-l) 



3 — 3 

the Lemma [I] that to — e < aj™ 0+Ar , b( no+JV ^ < vn + 2s. Similarly, when 

b {no+N) = m („ +W) ; we Qbtain m _ £ < a {no+N-X) ^{no+N-X) < m + 2e . By 

repeated use of the Lemma [I] in a similar manner, we hnally obtain 

v z a^ o) < to + (2 W + 2 - 2)e < M < M 

v z 6,[™ o) < to + (2 W+2 - 2)e < M < M in °\ 

which contradicts the definition of M^ n °\ Thus, we have proved to = M. Hence 
all the values converge to the same value C 1 ' 2N . 

In the preface of this article, we pointed out that the discrete model has to 
maintain the mathematical structures of a continuous model in the process of 
ultradiscretization. When we take ultradiscrete limit of C, it is also a conserved 
quantity of the pBBS. In fact, 



JY 



,(0)1,(0) UD 



i=i 



N 

EH 0) +sf} 



(15) 



gives the double number of balls in the pBBS. The number of balls is, to be 
sure, a conserved quantity of the pBBS. 

We can construct other conserved quantities of the recurrence formulae (^|) 
by means of another inverse ultradiscretization. 

From (H), we have 



J^n+l) ^ g-1 £(n+l) = j\(n) /\ 
j\{n+l) y g-1 £(n+l) _ j^(n) v _g(n) 



(16) 



When we consider the inverse ultradiscretization of the above equation, we have 



(n+l), (n+l) 
" : °i+l 



(n) , (n) 

a} o) 



n 



(n+l) 



l,(n+l) 
i+1 



,(«) , ,(n) 



(17) 
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Thus, for arbitrary A, (A + a\ n) )(X + b[ n) ) = (A + a|" +1) )(A + &££ X) ) and wc 
find that 



N 



C n (X) :=l[(\ + 4 n >)(\ + b. 



(18) 



does not depend on n, which means that any symmetric polynomial with respect 
to {a l -™' ) } and {b^} does not depend on n. Therefore, the ultradiscrete limit of 
such symmetric polynomials gives 27V conserved quantities S% , S2 , ■ ■ ■ , S2N of 



If we denote B 



(n) _ An) 
i N+i 



these conserved quantities are explicitly given 



a.s 



51 :— max 

i 

52 '■= max 

i<j 



in) 



A 



2N 



~>2N 



J2 A 

i=l 



(n) 



4 pBBS and combinatorial R matrix 

4.1 pBBS as periodic crystal lattice 

The BBS (of infinite number of boxes) has recently been reformulated from the 
theory of crystal and the combinatorial R matrix|| O, In this approach, a 
time evolution pattern of the BBS corresponds to a ground state configuration 

of a solvable lattice which has a symmetry of quantum algebra U' q {A^) 
The Boltzmann weight on every vertex of the lattice is given by combinatorial- R 
matrix of U' q (A^ ), and the states on each link are represented as the M-fold 
symmetric tensor representation Bm- For the simplest BBS in which only one 
kind of balls exist and all the box capacity is one, the lattice model has the 
space of horizontal links B1® 00 = ■ ■ ■ B± (g) B% <8> • • • , and that of the vertical 
links Boo® 00 = • • • (g) Boo ® Boo ® ■■ ■ ■ Here B M is understood as Boo = Bn 
(N ^> 1). Precisely speaking, N can be any positive integer which is greater than 

the number of the balls in the BBS. The combinatorial R matrix of U 1 q {A^) 
gives an isomorphism Boo <S> B\ — > B\ <£> Boo- The initial condition of the BBS 
corresponds to the initial state of the horizontal links of the lattice model. Since 
the number of balls is finite, the initial state •■•(£> \u 1 ^^) ® l^^ ) ® <8> • • • 

(G Bi®°° = ■ ■ ■ ® B x ® B x ® ■ ■ ■ ) satisfies the condition: 

K=°) = |0) for \n\ » 1, 

where |0) denotes the highest weight vector of B\. (The basis of B\ will be 
denoted by |0) and |1), and |0) corresponds to an empty box and |1) corresponds 
to a box with a ball.) The boundary condition for the vertical links is expressed 
as = |{0}) for (|n| ^> 1), where |{0}) is the highest weight vector of Boo. 
In general, we can replace Bi with Bg n and Boo with B Kt , where 9 n is the 
box capacity of nth box, and K t is the carrier capacity of tth carrying cart|l2|. 
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The solitonic natures of the BBS can be proved algebraically with the above 
settmg|l3| and the BBS can be extended to other quantum algebras [[u]. 

The pBBS discussed in the previous sections is also reformulated as a com- 
binatorial R matrix lattice model with periodic boundary condition. For the 
original BBS, time evolution is given by the isomorphism: 

T: Boo <8> Bf N - Bf N <g> Boo 

T: |{0})®|c(t))^| c (t + l))®|{0}) 

where \c(t)) £ Bf N is the state corresponding to the BBS at time t. For 
the pBBS, we have to take the trace of the vertical state, i.e., by regarding 
T e End EndB ®NBoo, we define the matrix T := Tr Boo T E Endc-Bf N , which 
gives a time evolution as: 

T:\c(t))~\c(t+1)). 

At a glance, one may think that \c(t + 1)) becomes a linear combination of 
many of the tensor products of B\ crystals, however, one tensor products of B\ 
crystals map to a unique tensor product of B\ and the resultant state exactly 
corresponds to the state of the pBBS at time step t+1. Even for the pBBS with 
M (M > 2) kinds of balls and various box capacities, the above lattice model is 
also well defined as far as the dimension of the vertical crystal is large enough, 
that is, Kt S> 1 for the vertical crystal B Kt of U' q (A^_ 1 ). We will show a proof 
of this fact for the case with one kind of balls. Since the evolution rule for M 
kinds of balls is decomposed into M steps as far as n t >• 1 , and only one kind 
of balls are moved at each step according to the same evolution rule, the proof 
is also true for the case with many kinds of balls. When k is small, however, 
the above construction will not give a unique tensor product and will not define 
an evolution rule of the pBBS. 

Since we treat only one kind of balls, the states are represented by U' q {A^p) 
crystal. First we consider the isomorphism B K ® Be — Bg ® B K given by the 
combinatorial R matrix. A state b in B K is usually denoted by a single raw 
semistandard Young tableaux of length n on letters 1 and 2. Instead we denote 
b = (y,K — y) where y is the number of 1 in the Young tableaux. For (y, k — 
y) (g) (x, 6 — x) ~ (x', 8 — x') ® (y', n — y'), we have the relation []i"T| 

x = y — min[/«, x + y] + miii[0, x + y] (19) 
y = x + min[«, x + y] — min[#, x + y\. (20) 

For k > 9, the relation is explicitly written as 

(x + y < o) 

6 {6<x + y<n) (21) 
9 (k < x + y) 

x' = x + y-y' (22) 

Now let n (n = 1, 2, . . . , N) be the capacity of nth box, and Kt be the capacity 
of the carrying cart at time step t. The state at time step t is given by |e(t)) € 
Bg x (g) Bg 2 ® ■ ■ ■ ® Bg N . Since Bg n is a U' q (A^) crystal, a vector b n S Bg n is 




13 




represented as b n — (x n , O n — x n ), where x n corresponds to the number of the 
balls in the nth box. We denote a state b\ <g> 62 <8 • • • ® by [x\, X2, ■ ■ ■ ,xn] 
for bi = (xj, 9i — Xi) (i — 1,2, . . . , N). The combinatorial R matrix of U' q (A^) 
gives the isomorphism T: 

T : B Kt ® (B 8l <g> Bg 2 <g> ■ ■ ■ <g> Bg N ) ~ $ Bg 2 $ • • • <8> B 9jv ) $ B Kt 
([2/0] <8> [zi,^, ... ,a;^] ~ [a;i,a; 2 . ■ ■ ■ »a^v] ® [y D 
From Eq. (5^), we have the following recurrence equations: 
Vn = F(y n -i;x n ,G n ) 

(X n + Vn-l < On) 
(0 n < X n + Vn-l < &t) 
[Kt < X n + Vn-l) 

(x n + Vn-l < n ) 

(0 n <x n + y n _ a < n t ) (23) 

(« t < X n + Vn-l) 



y'o = 

We see that the function F(y; x, 0) is a piecewise linear and monotonically in- 
creasing function of y which satisfies F(y + l;x,0) — F(y;x,0) = or 1 and 
< F(y;x,0) < nt- Since y' is a function of yo, we denote it by y' — 
Fn(vo] {xi}, {Oi}) := (F o F o ■ ■ ■ o F)(y )- The function Fjy is also monotoni- 

s / 

AT times 

cally increasing piecewise linear function, and F]y(yo+l; {^i}; Fn(vo) { x i}, 

{Oi}) — or 1 and < Fjy < K t . Thus, for < yo < K t , there is one and only 
one integer y* or one and only one finite interval [y*, y*] (y*, y* € Z) where the 
identity yo = -Fjv(yo; {Oi}) holds. Furthermore, from Eqs. (|23]), we find 
{x^} do not vary for y* < yo < y* ■ Therefore we conclude that, for given {x{\ 
and {Oi}, there is at least one yo (0 < yo < «t) at which y is equal to yo and that 
{x'j} are uniquely determined and have the same values for yo which satisfies 
y' = yo- The above conclusion means that T := Tre T S End-B^ ® • • • ® -B^ 
maps a state b\ ® • • • ® 6jv to the state which is also described by one tensor 
product of crystal basis. We summarize the above statement as a Theorem. 

Theorem 2 

If k is greater than any Oi (i — 1, 2, . . . , iV), the map T := Trg^T € EndBgj <8 
Bg 2 ®- ■ -®Bg N sends a tensor product of crystal basis to a unique tensor product 
of crystal basis of U' q (A^). Furthermore, for sufficiently large k, the statement 
holds for U'q(A^j) with arbitrary positive integer M. 

From the construction of the map, it is clear that the pBBS discussed in the 
previous section corresponds to the case Oi = 1 (i = 1,2,... , N), and we use 
this map to construct the pBBS with arbitrary box capacities and ball species. 

4.2 pBBS as A^_ x crystal chains 

When there are M balls of the same kind and N boxes, we can also reformulate 
the pBBS in terms of the combinatorial R matrix of U 1 q (Aw_i) and symmetric 
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A' 



tensor product Bm and Bm> where M' :— Oi — M (Fig. [|). 

i=l 

A crystal 6 G Bm can be denoted by b = (#1, £2, . . . , xn) with < Xi < 9i, 
^2iLi x i ~ M. We associate a state of pBBS with the crystal b in which Xi is 
the number of balls in the ith box of the state. For the crystal b, we define the 
dual crystal b — (x\,X2, ■ ■ ■ ,5jv) € Bm 1 , where Xi — 6i — x% (i = 1,2,--- ,N). 
Then the crystal 6' £ Bm associated with the state at time t + 1 is given by the 
combinatorial R matrix which gives the isomorphism Bm' ® Bm — Bm ® Bm> 
as 



R 



b®b^b' <g> &'. 



(24) 



From rcf. we see that this gives the same time evolution of the pBBS 

discussed above. As is shown in Fig. (j^), the time evolution is described in two 
twisted chains of Bm and Bm> ■ Note that by changing the crystal b and/or the 
b with another crystal (say B type crystal), we obtain other types of pBBS with 
a time evolution rule given by the isomorphism R. We may find interesting 
features in these CAs, however, the investigation of these pBBSs is a future 
problem. 



The isomorphism (24) has been shown to be expressed as an ultradiscrete KP 
equation (Eqs. (22) and (23) in ref. |jl2j), which serves another reason why we 
claim the pBBS is an integrable CA. Here we do not repeat the results in ref. ]l2]| , 
but show a similar formulae to (^). The space of the states, however, is no longer 
a finite field but 1i N . For X — (xi,X2, ■ ■ ■ , xn), Y — (yi, j/2, ■ ■ ■ , Vn) € Z , we 




Figure 5: Twisted chains of crystal A < ^ > _ 1 and pBBS 
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define max and min: Z x Z — * Z as 

(min[X,y]) i = minfajj.i/i] 
(max[X,Y]) i = max[xi,yi]. 

We also define the rotate shift to the right S: 1 N — > Z^ as 

SX = (xjv, 5Cl, X 2 , . . . jXiV-l). 

Let 9i(G Z >0 ) be the capacity of ith box and Xi(t) (0 < Xi(t) < 6i) be the 
number of balls in zth box at time step t. We denote the state of the pBBS at 
t by X(t) := (xi(t),X2(t), . . . , Xj^(t)). The state at t + 1, X(i + 1), is obtained 
from the following Theorem: 

Theorem 3 

Let A<® = X(t) and B^ = SX(t). We define and (n = 1, 2, . . . ) by 
the recurrence equations: 

f := min [A(") + B (n) , 0] , , 

\ B(" +1 ' :=5max[AW+B(")-0,O] ' ( 5j 

where 9 := (61,62, ■ ■ ■ ,6^). Then we obtain 

X{t + 1) = A {N -V - X{t), B {N -^ = 0. (26) 



Proof 3 

The proof of Theorem |^ is similar to that of the Theorem 1 and we simply show 
its outline. We define (i = 1,2,... , N) as in the proof of Theorem 1. Then 
we have the decomposition: 



N-l N-l 



X(t)=^D^, X(t + 1) = S l D (t) - (27) 

i=l i=l 

From the recurrence equations ( p^J and the properties of the similar diagram 
of evolution patterns, we can inductively show 

A (k) = A (k-i) + gk D {k) 

s -k-i B (k) = S -k B (k-i) _ D (k) ^ (28) 
gfc+ijjCfc+i) = min [ _ A (fe) t 

and we have 



AW = X(t) + ^5*I>W, (29) 

i=l 

S (n) = 5 n+l I ^ D (i) j 5 ( 30 ) 



Vz— n+1 



which completes the proof. 



16 



5 Summary 



In this paper, we introduced the pBBS which is an extension of the original 
BBS to that with a periodic boundary condition. We showed that the evolution 
rule of the pBBS is given by a recurrence Boolean formula which is regarded as 
an ultradiscretized algorithm of the calculation of 27Vth root of a given number. 
We also gave the conserved quantities of the recurrence formula. The relation 
to the combinatorial R matrix of U' q (A^_ 1 ) was clarified, and generalization 

of the pBBS with the symmetric tensor product representaion of U' q (A^_ 1 ) 
crystals was discussed. 

Since the pBBS takes only finite number of states, it has a fundamental 
cycle. Determination of this fundamental cycle is one of the future problems. 
In addition, an integrable equation usually has quasi-periodic solutions given by 
Theta functions. A state of the pBBS is expected to be obtained from the quasi- 
periodic solution through ultradiscretization. To obtain rigorous expression of 
the solutions to the pBBS through ultradiscretization is another future problem. 
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